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Abstract 

Low-energy electron impact excitations of S2 molecules are studied using the fixed- 
nuclei R-matrix method based on state-averaged complete active space SCF orbitals. 
Integral cross sections are calculated for elastic electron collision as well as impact 
excitation of the 7 lowest excited electronic states. Also, differential cross sections 
are obtained for elastic collision and excitation of the a^A^, and B^Ti~ states. 

The integrated cross section of optically allowed excitation of the B^Ti~ state agrees 
reasonably well with the available theoretical result. 
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1 Introduction 



Diatomic sulfur, S2, has been subject of many theoretical and spectroscopic 
investigations for long time [T|2|3] . We can find S2 molecules at various natural 
and industrial plasmas containing sulfur compounds. For example, emissions 
and absorptions of S2 molecules have been observed in the atmospheres of 
Jupiter [1] and its satellite lo [S]. They are also observed in the atmospheres 
of some comets [6]. In industrial condition, S2 molecules can be seen in re- 
active ion etching process using SFg molecules Sulfur lamps contain S2 
molecules as an important ingredient [8j. Although electron collision with S2 
molecules is an important elementary process in these plasmas, there has been 
little work on this subject. As far as we are aware, no experimental measure- 
ment nor theoretical calculation of electron-S2 elastic cross section have been 
performed. Garrett et al. ^ calculated integral cross section for electron im- 
pact excitation of the S2 B^E~, 2^E~ and B"^n„ states using the semiclassical 
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impact-parameter (IP) method extended to include nuclear motion. Since the 
IP method is designed to treat optically allowed transitions, electron impact 
excitation to the other electronic state was not investigated. Le Coat et al. 
[TU] performed experimental measurement on dissociative electron attachment 
of S2 molecules and identified two resonances, however, they did not present 
absolute value of the cross section. 

Recently, we have performed the R-matrix calculations on the spin-exchange 
effect in electron collision with homo-nuclear open-shell diatomic molecules 
including S2 [H]. We have also calculated the integral cross section of elastic 
electron collision with S2 molecules for the first time. Here, we extend our 
previous study to the electron impact excitations of S2 molecules. Since S2 
molecule has a valence electron structure similar to O2, comparison of the 
cross sections will be interesting. Also, the cross section of electron impact 
excitation of the S2 B^E~ would be important in analyzing sulfur plasma. As 
in our previous works on electron impact excitation of O2 and N2 p!^13fl4j . 
we employ the fixed-nuclei R-matrix method using state-averaged complete 
active space SCF molecular orbitals. The size of basis set for the scattering 
electron is shghtly extended in this work compared to our previous calculation 
on the spin-exchange effect [TT]. 



2 Theoretical methods 

The details of the R-matrix method has been described in the literature 
p!5|16|17j . thus we do not repeat general explanation of the method here. 
We used a modified version of the polyatomic programs in the UK molecular 
R-matrix codes [T3| in this work. These programs utilize the Gaussian type 
orbitals (GTO) to represent target electronic states as well as a scattering elec- 
tron. In the present R-matrix calculations, we have included 13 target states; 
X^E-, a'Ag, b'E+, c^S-, A'^A^, A^SJ, B'%, E^S", l^n,, l^A,, S"'n„, 
1^S+, and l^n^. These target states were represented by valence configura- 
tion interaction wave functions constructed by state averaged complete active 
space SCF (SA-CASSCF) orbitals. In this study, the SA-CASSCF orbitals 
were obtained by calculations with MOLPRO suites of programs [TS]. The 
target orbitals were constructed from the cc-pVTZ basis set[T9]. Fixed-nuclei 
approximation was employed with inter- nuclear distance of 3.7 ao. The radius 
of the R-matrix sphere was chosen to be 13 ao in our calculations. In order to 
represent the scattering electron, we included diffuse Gaussian functions up to 
/ = 5, with 13 functions for Z = 0, 11 functions for Z = 1, 10 functions for I = 
2, 8 functions for / = 3, 6 functions for / = 4 and 5 functions for / = 5. Expo- 
nents of these diffuse Gaussians were taken from Faure et al. [20]. In addition 
to these continuum orbitals, we included 8 extra virtual orbitals, one for each 
symmetry. The construction of the configuration state functions (CSFs) for 



2 



the electron- molecule system is the same as in our previous paper [TT]. Note 
that The R-matrix calculations were performed for all 8 irreducible represen- 
tations of the D2/1 symmetry, Ag, B2U, B^u, B^, Bi^, B^g, B2g and A^, in 
doublet and quartet spin multiplicities of the whole system. 

One of the transitions studied in this work, the excitation of the B^'E~ state 
from the ground state, is optically allowed. Thus, we have to consider the 
effect of transition dipole moment between these two states. A lot of / partial- 
waves has to be included in the R-matrix calculation to obtain converged 
cross sections because of the long-range interaction of the dipole, although it 
is difficult to include partial waves with / > 7 in the usual ab initio R-matrix 
calculation. In this work, the R-matrix calculations are performed with partial 
waves up to / = 5. The effects of the higher / partial waves are included by 
the Born closure approximation as in the previous works pill^ . Following 
Gibson et al.[22j, we evaluate the differential cross sections (DCSs) with the 
Born correction, da^^ /(El, by the expression, 



da^_da^ 



^^R-matrix ^^FBA 

In dn~ 



Here, da^^^ / dVL is the DCS obtained by the first Born approximation, da^~^^^^^^ / dVL 
is the cross section obtained by the R-matrix calculation and da^^ /dVL is the 
DCS from the first Born approximation including the same number of partial 
waves as in the R-matrix calculation. da^^/dVL is evaluated by the angular 
momentum representation of the T-matrix elements for the first Born approx- 
imation. These T-matrix elements as well as da^^^/dVL are available in close 
form [23]. The total cross sections are obtained by the integration of eq.(l). 



3 Results and Discussion 



In table [H excitation energies of S2 molecule obtained from the CASSCF cal- 
culation in this work are compared with MRD CI vertical excitation energies 
of Hess et al. [21], MRCI adiabatic excitation energies of Kiljunen et al. [25] 
and experimental values quoted in Hess et al. Our results agree well with the 
previous calculations and experimental results for the lowest two excitations. 
For excitations to the higher electronic states, deviations in excitation energies 
become larger partly because of difference of adiabatic and vertical excitation 
energy. 

In figure [H the integral cross sections (ICSs) are shown for electron S2 elastic 
collision and excitations of the a}/S.g, and c^E~ states. The magnitude 

of the elastic ICS is about 20 xlO^^^cm^ in low energy region below 3 eV, 
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then it increases to be 30 xlO~^^cm^ at energies above 5 eV. The ICS of 
the excitation to the a}^g state increases gradually from threshold to 8.5 eV, 
where it takes a maximum value of 0.35 xlO~^^cm^, then it decreases again. 
The ICSs of the excitation to the fe^S^ and c^S^ states also increase gradually 
from threshold. In both cases, the maximum value of the ICS is about 0.1 
xlO~^^cm^. Around 2.7 eV, a sharp resonance peak with width 0.08 eV exists 
in the ICSs of the and excitations. Also, a kink structure is observed 
in the elastic ICS at the same energy. All of these structure belongs to the 
symmetry partial cross sections. We analyzed the CSFs and found that the 
kink and peaks at 2.7 eV are likely related to a resonance with configuration 
(core)^°(4crg)^(4cr„)^(5crg)^(27r„)^(27rg)'^, which is obtained from an attachment 
of the scattering electron to the excited c^E~, A'^A^ and states of S2 

with configuration {coieY^{4:agf{AauY{bagf{2'Kuf{2TTgf. The width of this 
resonance is about 0.08 eV at R = 3.7 ao. We checked the behaviour of this 
resonance as a function of bond-length, and found that it approaches the 
atomic limit S('^P) +S~(^P). Since the location of this resonance is higher 
than the atomic limit, dissociative electron attachment may occur through 
this resonance. 

In figure [21 the ICSs for electron S2 collisions are shown for excitations to the 
A'^K, B'^Iig and B^T.- states. The ICS of Garrett et al. obtained 

by the impact-parameter method is also compared with our excitation ICS to 
the B^T,~ state. In general, the slopes of these ICSs near threshold are steeper 
than those in figiH The maximum values of the ICSs below 15 eV are about 
0.3, 0.1, 0.6 and 3.0 for the excitation to the A'^A„, ^^2+, B'^Iig and B^T.- 
states, respectively. For the optically allowed B^T,^ state excitation, the ef- 
fect of the high / partial waves is included by the Born closure approximation 
formula given in eq.(l). The magnitude of the Born correction is small below 
10 eV, however, it increases as the scattering energy increases. The fraction of 
the correction to the R-matrix cross section becomes about 20% at 15 eV. Our 
ICS with the Born correction and the previous result of Garrett et al. [9] agree 
reasonably well above 8 eV. We used the CASSCF value for the B^Y^~ excita- 
tion energy in this work, whereas Garrett et al. [9J employed the experimental 
value for the excitation energy. In contrast to the fixed-nuclei approximation 
in our calculation, Garrett et al. [9j included the effect of S2 vibration in their 
calculation. Because of these differences, the results of Garrett et al. [9] and 
our ICS do not agree well near the excitation threshold. A small peak is seen 
at 14.2 eV in the ICSs of the A'^A„ and ^^S^ state excitations, which is 
originated from the '^11^ symmetry partial cross sections. Since the location of 
this peak is higher than the highest energy S2 electronic state included in the 
present R-matrix calculation, we cannot determine whether this peak belongs 
to the real resonance or pseudo-resonance. 

In figur^ni the differential cross sections (DCSs) are shown for elastic electron 
S2 collisions as well as excitations to the a^A^,, and states. The 
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elastic DCSs are enhanced in forward direction and tend to be more forward- 
enhanced as the scattering energy increases. In contrast to the elastic DCSs, 
the excitation DCS to the a^A^ state has backward-enhanced character in 
general. However, the magnitude of forward scattering cross section increases 
as the scattering energy increases from 7 to 13 eV. Our excitation DCSs to the 
state approach zero near and 180 degrees, because of a selection rule 
associated with S^-E" transition [26|l27j . Around 90 degrees, the magnitude of 
the DCS is about 0.06~0.08 xlO~^^cm^sr~^ and does not depend much on the 
scattering energy. The B^T,^ state excitation DCSs are shown in fig. [3] panel 
(d). In addition to the R- matrix results, the DCSs with the Born correction 
obtained by eq.(l) are also shown for scattering angles below 25 degrees. For 
larger angles, the magnitude of the correction is expected to be small and not 
shown here. The contribution of the Born correction to the DCS is small for 
7 eV, however, it dominates the total DCS near zero degree at 10 and 13 eV. 
Although the magnitude of the R-matrix DCS at forward direction decreases 
as the scattering energy increases from 7 to 13 eV, the DCS with the Born 
correction at forward direction increases as energy increases. 

The elastic and excitation cross sections of the a^A^ and states in elec- 
tron S2 collisions are about two times larger than corresponding cross sections 
in electron O2 collisions studied in our previous paper [12]. Although the ^Uu 
resonance peaks can be observed in the a^Ag and excitation cross sec- 
tions in both e-02 and e-S2 collisions, the width of the peak is much broader 
in e-02 case. In e-02 elastic collision, a narrow ^II^ resonance peak is seen 
below 1 eV. In e-S2 elastic case, the energy of the S2 state is stabihzed 
below the energy of the S2 ground state and cannot be observed in the cross 
section. Other than these resonance features, the profiles of the cross sections 
are similar in e-02 and e-S2 collisions. 

In this work, we employed the fix-bond approximation for the R-matrix cal- 
culation. In our previous studies on electron impact excitations of O2 and N2 
molecules |12lll3lll4j . we also used the same fix-bond approximation and got 
good agreement with available experimental results, even if the positions of 
the potential curve minimum are different between the ground state and the 
excited state. Thus, the results of this study is also expected to be accurate 
enough. For more precise comparison of the excitation cross section of the 
i?^E~ state with the previous results of Garrett et al. inclusion of the 
vibrational effect may be necessary. Such kind of calculation is possible by the 
non-adiabatic R-matrix method or the adiabatic averaging of the T-matrix 
elements, and will be performed in future. 

We have carried out the R-matrix calculations with the maximum / quantum 
number 4, 5 and 6 to check convergence. Except for the -B^S~ state excitation, 
the ICSs and DCSs are converged at 1=5. For the B^T,' excitation, the ICS 
with the Born correction is also converged at 1=5. For the B^H' excitation 
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DCSs with the Bom correction, however, the convergence is achieved only 
below 25 degrees. Note that the similar situation was observed in Gibson et 
al.[22]. Although the S^S" excitation DCS is not converged above 30 degrees, 
the effect of the higher / partial waves is expected to be small because the 
magnitude of the Born DCS itself is small at larger scattering angles. In prin- 
ciple, this convergence problem can be solved [28H29] by applying the Born 
correction to T-matrix elements, where we applied the correction at the R- 
matrix DCS in this work. We will calculate this kind of Born correction at 
T-matrix level when more accurate excitation DCS is required in future. 



4 Summary 

In this work, we have studied the low-energy electron impact excitations of 
S2 molecules using the fixed-bond R-matrix method based on state-averaged 
CASSCF molecular orbitals. Thirteen target electronic states of S2 are in- 
cluded in the model within a valence configuration interaction representations 
of the target states. Integral cross sections are calculated for elastic electron 
collision as well as impact excitation of the 7 lowest electronic states. Also, 
differential cross sections are shown for elastic collision and excitation of the 
a^Ag, b^^^ and B^J]^ states. For the excitations of the a^Ag and states, 
a narrow resonance peak is observed in the ICSs at 2.7 eV. For the elas- 
tic and the excitation collisions of the a^Ag and states, the shapes of 
the cross sections are similar to those in electron O2 collisions, however, the 
magnitudes of the cross sections are two time larger in electron S2 collisions. 
Our ICS of the B^T,~ state excitation agrees reasonably well with the previous 
result of Garrett et al. [S]. However, the ICS near threshold does not agree 
well, because of difference in excitation energy employed in calculation as well 
as the treatment of the vibrational effect. 
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Table 1 

The vertical excitation energies of the first 8 excited states for S2 molecule, with 
the previous MRD CI results of Hess et al.[23], MRCI results of Kiljunen et al.[25] 
and experimental values quoted in Hess et al.[2l]. The unit of energy is eV. 



State 


This work 


Previous MRD CI 


Previous MRCI 


Expt. 




0.00 


0.00 


0.00 


0.00 


a^Ag 


0.60 


0.68 


0.55 


0.71 




0.92 


1.04 


0.99 


0.99 




2.77 




2.45 






2.93 




2.59 






3.03 




2.58 




B''% 


4.84 


4.63 


4.36 


4.38 




5.03 




3.89 
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Electron impact energy (eV) 

Fig. 1. Integral elastic (panel a) and excitation cross sections of the a^Ag (panel b), 
(panel c) and c^S~ (panel d) states. 
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Electron impact energy (eV) 

Fig. 2. Integral excitation cross sections of the ^''^A^ (panel a), A'^T,^ (panel b), 
B''^Ilg (panel c) and B^T,~ (panel d) states. The R-matrix cross section with Born 
correction and the previous results of Garrett et al. [9] are also included in panel d. 
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90 120 150 180 



Fig. 3. Differential cross sections of elastic collision (panel a) and excitation of the 
a^Ag (panel b), (panel c) and B'^T,~ (panel d) states. The R-matrix cross 

sections with Born correction are also included in panel d. 
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